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1. Introduction 

The production of heavy quarks in high energy 
processes has become an increasingly important 
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subject of study both theoretically and experi- 
mentally. The theory of heavy quark produc- 
tion in perturbative Quantum Chromodynam- 
ics (PQCD) is more challenging than that of 
light parton (jet) production because of the new 
physics issues brought about by the additional 
heavy quark mass scale. The correct theory must 
properly take into account the changing role of 
the heavy quark over the full kinematic range 
of the relevant process from the threshold region 
(where the quark behaves like a typical "heavy 
particle") to the asymptotic region (where the 
same quark behaves effectively like a parton, sim- 
ilar to the well known light quarks {u, d, s}). 

We review theoretical methods which have 
been advanced to improve existing QCD calcula- 
tions of heavy quark production, and the impact 
on recent experimental results from HERA and 
the Tevatron. 

The ACOT renormalization scheme provides a 
mechanism to incorporate the heavy quark mass 
into the theoretical calculation of heavy quark 
production both kinematically and dynamically. 
In 1998 Collins^ extended the factorization the- 
orem to address the case of heavy quarks; this 
work provided the theoretical foundation that al- 
lows us to reliably compute heavy quark processes 
throughout the full kinematic realm. 
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Figure 1. Characteristic Feynman graphs which 
contribute to DIS heavy quark production: a) the 
LO C(ag) quark-boson scattering QV — > Q, 
b) the NLO 0{ag) gluon-boson scattering gV 
QQ, and c) the NNLO C(q;|) boson-gluon scat- 
tering gV -> gQQ. 

1.1. NLO DIS calculation 

Figure [l] displays characteristic Feynman 
graphics for the first few orders of DIS heavy 
quark production. If we consider the DIS produc- 
tion of heavy quarks at 0{ag) this involves the 
LO QV -> Q process and the NLO gV QQ 
process. 

The key ingredient provided by the ACOT 
scheme is the subtraction term (SUB) which re- 
moves the "double counting" arising from the re- 
gions of phase space where the LO and NLO con- 
tributions overlap. Specifically, the subtraction 
term is: 



CrSUB = /o «> Pg 



a SUB represents a gluon emitted from a pro- 
ton (/g) which undergoes a collinear split- 
ting to a heavy quark (Pg^g) convoluted with 
the LO quark-boson scattering (Jqv^q- Here, 
Pg^gi^^fJ-) = §i Pg^cjx) where 

Pg-,c{x) is the usual MS splitting kernel. 

1.2. When do we need Heavy Quark PDFs 

The novel ingredient in the above calculation 
is the inclusion of the heavy quark PDF con- 
tribution which resums logs of ln(/i^/mQ). One 
can ask the question: When do we need to con- 
sider such terms? The answer is illustrated in 
Figure [2] where we compare the DGLAP evolved 
PDF fc{x,iJ,) with the single splitting perturba- 
tive result 
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Figure 2. Comparison of the DGLAP evolved 
charm PDF fc{x, fj,) with the perturbatively com- 
puted single splitting (SUB) fc{x, fj.) = fg{x, /x) (g) 
Pg~^c charm evolution vs. ^ in GeV for two rep- 
resentative values of x. 



The DGLAP PDF evolution sums a non- 
perturbative infinite tower of logs while the 
SUB contribution removes the perturbative single 
splitting component which is already included in 
the NLO contribution. Hence, at the PDF level 
the difference between the heavy quark DGLAP 
evolved PDF fq and the single-splitting pertur- 
bative /q will indicate the contribution of the 
higher order logs which are resummed into the 
heavy quark PDF. Here, we shall find it conve- 
nient to define fq = fg® Pg^Q which represents 
the PDF of a heavy quark Q generated from a 
single perturbative splitting. 

For jjL ~ mg we see that fq and fq match quite 
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closely, whereas /q and /g differ significantly for 
IJ, values a few times mg. While the details will 
depend on the specific process, in general we find 
that for fj, scales 3 to 5 times mq the terms re- 
summed by the heavy quark PDF can be signifi- 
cant. 

2. The ACOT Renormalization Scheme 

2.1. Massive vs. Massless Evolution 

Another useful result that arises from the proof 
of Collinsjl] is that we can use mass-independent 
(massless) evolution kernels to evolve the heavy 
quark PDFs without any loss of accuracy as com- 
pared to a mass-dependent (massive) evolution 
kernel. [2] Specifically, Collins demonstrated that 
consistent application of the formalism correctly 
resums the massive contributions up to higher- 
twist corrections ©(Ag^^/Q^) and that there are 
no errors of order ©(mg/Q^). 

This result is illustrated in Figure |3] where we 
compare the results of a NLO DIS heavy quark 
production calculation using massless and mas- 
sive DGLAP evolution kernels. In Fig. |3^) we 
see that while the choice of massive or massless 
kernels significantly changes the individual LO 
and SUB contributions, the difference LO — SUB 
which contributes to the total {TOT = LO — 
SUB + NLO) is minimal. This numerically ver- 
ifies that the choice of massive or massless evolu- 
tion kernels is purely a scheme choice which has 
no physical content. 

While we see this result demonstrated numer- 
ically in Figure |3j the underlying reason for this 
result is closely related to the previous observa- 
tions made regarding Figure [2] The LO result is 
given by LO ~ /g (8i ctq^q and the subtraction 
term is given by SUB ~ /g (8) Pg^Q ® uq^q. If 
we expand the DGLAP equation for /g in the re- 
gion /I - mg we find fg y fg Pg^g + C(a^); 
thus, we have LO ~ /g (g) Pg_^g ® crg_+g -|-0(a|). 
We observe that while LO and SU B individually 
depend on the specific splitting kernels, the com- 
bination LO — SUB is insensitive to whether we 
use the massive or massless kerneljf] 

^While we have given a heuristic description of this re- 



Therefore, we conclude that so long as the 
splitting kernels Pa~^b are matched between the 
DGLAP evolution and the definition of the sub- 
tractions (SUB), the choice of a massive or mass- 
less DGLAP evolution kernel was purely a choice 
of scheme and the physical results are invariant. 

2.2. S-ACOT 

In a complementary application, it was ob- 
served that the heavy quark mass could be set 
to zero in certain pieces of the hard scattering 
terms without any loss of accuracy. This modi- 
fication of the ACOT scheme goes by the name 
Simplified- ACOT (S-ACOT) and can be summa- 
rized as follows. 

S-ACOT: For hard-scattering processes with in- 
coming heavy quarks or with internal on- 
shell cuts on a heavy quark line, the heavy 
quark mass can be set to zero (mg = 0) for 
these pieces. |3] 

If we consider the case of NLO DIS heavy quark 
production, this means we can set mg = for the 
LO terms {QV Q) as this involves an incoming 
heavy quark, and we can set rag — for the SUB 
terms as this has an on-shell cut on an internal 
heavy quark line. Hence, the only contribution 
which requires calculation with mg retained is 
the NLO gV QQ process. 

Figure |4] displays a comparison of a calculation 
using the ACOT scheme with all masses retained 
vs. the S-ACOT scheme; as promised, these two 
results match throughout the full kinematic re- 
gion. 

2.3. ACOT-x 

In the conventional implementation of the 
heavy quark PDFs, we must "rescale" the Bjorken 
X variable as we have a massive parton in the fi- 
nal state. The original rescaling procedure is to 
make the substitution x ^ x{l + m^/Q^) which 
provides a kinematic penalty for producing the 
heavy charm quark in the final state. [3] As the 
charm is pair-produced by the g ^ cc process, 

suit (in which we used some illustrative approximations), 
we emphasize the proof applies to all cases and does not 
require any such approximations. 
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Figure 3. Comparison of heavy quark DIS structi 
independent (massless) evolution. 
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Figure 4. Comparison of schemes for NLO DIS 
heavy quark production as a function of Q. We 
display calculations using the ACOT, S-ACOT, 
Fixed-Flavor Number (FFN), and Zero-Mass 
Variable Flavor Number (ZM-VFN) schemes. 
The ACOT and S-ACOT results are virtually 
identical. 

there are actually two charm quarks in the final 
state — one which is observed in the semi-leptonic 
decay, and one which goes down the beam pipe 
with the proton remnants. Thus, the appropri- 
ate rescaling is not x x{l + m^/Q^) but in- 
stead X ^ X — + (2mc)^/(5^); this rescal- 
ing is implemented in the ACOT-x scheme, for 
example. |'5i6l7] The factor (1 -I- {2mc)^/Q^) rep- 
resents a kinematic suppression factor which will 
suppress the charm process relative to the lighter 
quarks. 
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Figure 5. Calculation of DIS heavy quark pro- 
duction for a variety of schemes. 

2.4. Numerical Comparison 

Having introduced the various theoretical is- 
sues which enter the calculation of the heavy 
quark process, we illustrate the numerical size of 
these choices for the case of DIS heavy quark pro- 
duction. 

In Figure |5] we display the charm structure 
function F2{x,ii) for a variety of schemes and 
orders. LO represents the 0(a°) QV Q pro- 
cess. NLO includes the 0{al) processes (primar- 
ily gV — *■ QQ) in the massless approximation. In 
the Fixed-Flavor-Scheme (FFS) the heavy quark 
PDF is set to zero; hence, at 0{al) this only re- 
ceives contributions from gV QQ. The ACOT 
and S-ACOT schemes are virtually identical — the 
curves are indistinguishable in this plot. Finally, 
the implementation of the ^-prescription for the 
S-ACOT scheme (the ACOT-x would yield iden- 
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Set 


# points 


CTEQ6HQ 


CTEQ6M 


6M®GM 


6HQ«)ZM 


ZEUS 


104 


0.91 


0.98 


2.84 


3.72 


HI 


484 


1.02 


1.04 


1.50 


1.22 


TOTAL 


1925 


1.04 


1.06 


1.26 


1.30 



Table 1 

Table of per point for the individual HERA data sets, and for the TOTAL of all data sets. (Non-HERA 
data sets are not displayed.) The results are shown for CTEQ6HQ PDF using the General Mass (GM) 
ACOT scheme, and CTEQ6M PDF using the zero-mass (ZM) MS scheme. We note the increased for 
mixed schemes using CTEQ6M with the GM ACOT scheme, and the CTEQ6HQ with the ZM scheme. 



tical results) provides some additional suppres- 
sion in the region fj, ~ mq. To this order, our 
best theoretical estimate of the true cross section 
would be either the ACOT-x or equivalently S- 
ACOT-x. 

To see the effect of these different results in the 
context of a global fit we display the results for the 
CTEQ6M and CTEQ6HQ PDFs sets in Table [l] 
Both the fits using a consistent application of the 
ACOT and MS schemes yield good results. In 
contrast, if we mismatch the scheme used in the 
PDF with that in the cross section calculation 
we observe a dramatic increase in the values 
obtained. This result underscores the importance 
of using properly matched calculations. 

2.5. Heavy Quarks at the Tevatron 

In the previous discussion we have primarily 
focused on DIS production of heavy quarks for 
illustrative purposes as the formalism is easier to 
layout when there is only a single hadron in the 
initial state. Nevertheless, the same principles 
that we have used in the the DIS case can be 
applied to that of the hadron-hadron initial state 
as appropriate for the Tevatron and the LHC. 

Historically, the predictions of b-production 
at hadron-hadron colliders have been a chal- 
lenge; the early results from the Tevatron were 
a factor of 2 to 3 larger than the theoretical 
predictions. NLO QCD corrections to the LO 
99 ~^ QQ process were formidable and yielded 
large corrections. |11I12I13I14I15] It is interesting 
to observe that if the heavy quark PDF is taken 
into account so that the LO contribution con- 
sists of both gg QQ and gQ — > gQ, then 
the computed NLO contributions (with appro- 



priate subtractions) are thereby reduced suggest- 
ing improved convergence of the perturbation 
theory. |16I17|18. 19 20 .21] 

Ref. |8] performs a systematic comparison of 
the GM-VFNS and ZM-VFNS using results of 
an updated analysis of hadronic b-production at 
the Tevatron. Figure [6] displays these results for 
the Tevatron in the central rapidity region as 
compared with the CDF data.[9 lOj The result 
is that the finite mass effects moderately enhance 
the pt distribution in the region pr ~ 2m// by 
about 20%, and this enhancement decreases at 
larger pT- For intermediate to large pT values 
{pt > niH) the three calculations (GM-VFNS, 
ZM-VFNS, FFN) match quite closely, and are in 
good agreement with the data. Conversely, if we 
use the FFN result with the historic values for the 
PDF and as we find this prediction is roughly a 
factor of 3 below the data. 

The excellent agreement between data and the- 
ory for this process is an important achievement 
and represents the culmination of many years of 
effort by both the theoretical and experimental 
community. 

3. Schemes used for Global Analysis 

The ACOT scheme and variants were used 
for the CTEQ series of global PDF fitsQ For 
the MRST/MSTW series of global PDF fits the 
Thorne-Roberts (TR) scheme was used. As these 
two sets of PDFs are widely used it is of interest 
to compare and contrast these approaches. Fig- 
ure [7] displays a diagrammatic comparison of the 

■^Specifically, ACOT was used for CTEQ6HQ, and S- 
ACOT-x was used in CTEQ6.5 and CTEQ6.6. 
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Figure 7. Diagrammatic comparison of TR and ACOT type schemes for the case of DIS. This diagram is 
schematic to emphasize the similarities and differences. The leading-order (LO) process is a 0{a'^) boson 
scattering from a heavy quark, e.g. jQ — > Q; the NLO 0{al;) correction arises from — > QQ, and the 
NNLO 0{ag) correction arises from 717 QQg- 



TR |22I23] and ACOT type schemes. While these 
schemes may appear quite different at first glance, 
they differ by higher-order terms which will be re- 
duced as we increase the order of our perturbation 
theory. 

In perturbation theory, we compute our observ- 
ables to a fixed order N in as', hence, we trun- 
cate the perturbation expansion at 0{ag), and 
we have neglected terms of order 0{ag^^). In 
brief, the difference between these two approaches 
amounts to adding different ©(a^"*"^) higher or- 
der terms. Thus, these two approaches will agree 
on the contributions up to 0{ag). We will now 
review the motivation and consequences of adding 
the differing higher order terms. 



3.1. Leading-Order (LO) {aV) 

If we work at Leading-OrdeiP] (LO) ag, when 
the heavy quark PDF is an "active" parton (typi- 
cally /i > mn) the LO contribution is 'y + Q Q. 
However, when the heavy quark PDF is not 
an "active" parton (typically /i < mn) the LO 
contribution vanishes. For the ACOT scheme, 
no higher order terms are added to this results. 
Hence for scales /i < m//, the LO answer is 
zero and we expect large corrections to this re- 
sult at NLO. For the TR scheme, a portion of the 
75 ~^ QQ contribution is added; for /i < mn the 
full 75 QQ term is included, and for /x > itih 
the 717 ~* QQ term frozen at // = Q to avoid 
any difficulty with large logarithms of the form 
ln(m_f//^). 



^Here, we define the order of the calculation according to 
the power of as; thus LO is a^, NLO is a], etc. 
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Figure 6. From Ref. |8], the transverse mo- 
mentum distribution da/dpT for pp —^ BX at 
i/s=1.96TeV. The results are shown for the Gen- 
eral Mass (GM) Variable Flavor Number (VFN) 
scheme and the Zero Mass (ZM) Variable Fla- 
vor Number (VFN) scheme. Additionally, results 
are shown for the Fixed Flavor Number (FFN) 
scheme with both recent PDFs (dot-dashed line) 
and the historical PDFfe (dotted line). The data 
is from the CDF collaboration. |9I1Q] 

Consequently, in the fi < niH region the TR 
scheme yields a finite LO result while the ACOT 
scheme yields zero. While both schemes formally 
agree at 0{a^), clearly the 0{ag) terms can be 
important, particularly in the fj, < mu region. 

3.2. Next-to-Leading-Order (NLO) {a\) 

If we work at NLO (a^), for the low region we 
now include 75 QQ as well as the 7 + Q ^ Q 
process]^ If we again look in the region /i < niH, 

®Note, in Figure [t] and in the discussion the diagrams 
an processes are schematic and illustrative. For exam- 



we find that while the ACOT scheme yielded zero 
at LO, it now obtains a finite result at NLO. For 
the TR scheme, in addition to the above terms, a 
portion of the 73 gQQ contribution is added; 
again, for fi > niH the jg — > gQQ term is frozen 
at /i = Q to avoid any difficulty with large loga- 
rithms. 

As before, both the TR scheme and ACOT 
scheme formally agree at 0{ag), but they will dif- 
fer by the separate NNLO ^^(q;!) terms that have 
been included. In contrast to the LO case where 
the ACOT scheme yielded zero for /z < tuh, both 
schemes give finite results in all kinematic region; 
hence, the relative difference will be reduced. 

3.3. General Comparisons at Order 

Let us make some observations regarding these 
schemes at a general order in perturbation 
0{ag). We observe that for a given set of pro- 
cesses calculated to , we can implement the 
TR scheme to ©(a^"^) and the ACOT scheme 
to 0{ag). For example, at NLO we note that 
the ACOT scheme involves only graphs of or- 
der al while TR utilizes graphs of order a^. At 
present we know the 0(a|) massive neutral cur- 
rent process (717 gQQ, jQ ggQ and asso- 
ciated graphs); hence, this allows us to compute 
the TR scheme to 0{ag) and the ACOT scheme 
to ©(al). In contrast, the massive charged cur- 
rent process is known only to 0{a\;); hence, this 
allows us to compute the TR scheme to 0{a'g) 
and the ACOT scheme to 0{al). 

We note that recent improvements of theoret- 
ical techniques have enabled significant advances 
in the calculation higher-order heavy quark pro- 
cesses. For example, Ref. |25] has obtained the 
asymptotic results for F^'^{x, at the 3-loop or- 
der, and recently Ref. [26] has extended this work 
for the case of F^'^{x, fi). 

In general, the TR scheme achieves in practice 
the same highest asymptotic order as ACOT by 
some modeling of terms below Q^ — mg which 
become (relatively) unimportant at high Q^. As 

pie, at NLO we include both 7Q Qg and -yg QQ as 
well as all the corresponding subtractions. For details see 
Refs. [22 24 . 
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Figure 8. The b-quark PDF xfb{x,Q) with 
NNLO matching conditions for 3 choices of x. 

we move to higher order calculations, the differ- 
ences between these schemes will be reduced as 
they arise from uncalculated higher-order contri- 
butions. 

4. NNLO and Beyond 

Although NLO is the state-of-the-art for many 
calculations, improved experimental precision de- 
mands that we strive toward a NNLO accuracy. 
When we consider PDFfe for heavy quarks at 
NNLO, there are a number of new elements that 
enter. 

One consequence is that the PDFs are no 
longer continuous across the heavy flavor thresh- 
old. Even more, when matching charm and bot- 
tom across their thresholds, they start from nega- 
tive values as illustrated in Figure |8] The match- 
ing conditions have been computed by a number 
of groups. j27l28) and at NNLO PDFs will have 
discontinuities of order 0{ag) when we transi- 
tion from Np to Np + 1 flavors. While we may be 
uncomfortable with discontinuities in our PDFs, 
we are reminded that the PDFs are not physi- 
cal observables, but instead are only theoretical 
constructs which depend on (arbitrary) renormal- 
ization schemes and scales^ 

At NLO, the point /i — uiq is special because 

^Recall asifJ.) is also an unphysical theoretical construct; 

this has discontinuities across flavor-thresholds at order 
^3 



O" 











N,=3 






J 




Np=4 / 


N=5 1^ 


N=6 



Qn 



Q„ 



Np=3 1 1 



transition 
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Figure 9. The upper figure schematically rep- 
resents how each calculation with a set num- 
ber of flavors Np has a region of applicability. 
The transition from the Np-\ scheme to the Np 
scheme should be in the vicinity of the mTv^ mass, 
but need not occur exactly at = Trij^p. The 
lower figure illustrates that multiple PDFs can co- 
exist for ^ > rriMp with matching performed at 
H = rriNp . 

/^''{xjItiq) = mg); this is because the 

constant term in the matching equation happens 
to be zero at NLO. Because of this "accident" it 
was common to use ^ = mq as both the Matching 
Point and the Transition Point. 

At NNLO the point /i — mq no longer has 
these special properties as the transition from 
Np to Np + 1 will necessarily have discontinu- 
ities at any value of /i; hence, it may be desir- 
able to choose the Matching Point and the Tran- 
sition Point at different values of /i. As these two 
point are not usually distinguished, let us high- 
light their key features. 

Matching Point ^m- The value of where the 
Np -t- 1 scheme is defined in terms of the 
Np scheme by a relation of the form: 

f:^-+\x,^J) = A,,®f^-{x,^i). 

Transition Point ^.t- The value of where the 
user chooses to transition from the Np 
scheme to the Np -\- 1 scheme. 
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Figure [9] schematically represents how each cal- 
culation with a set number of flavors Np has a 
particular region of applicability where it is best 
suited to describe the "true" physics. The com- 
plete description of the physics throughout the 
full kinematic range will therefore consist of a 
patchwork of schemes which are "sewn together." 

The Transition Point: It is easy to imagine 
situations where we would not want to automati- 
cally transition between schemes at fi — mq. For 
example, consider we are analyzing data in the 
range /i S [2, 5] GeV. The bulk of the range is in 
the Np — 4 flavor region as > nric ^ 1.3 GeV, 
but a small portion of the range extends above 
the Np = 5 flavor region as /j, > rub ~ 4.5 GeV. 
In the region ^ G [4.5, 5] GeV it would be incon- 
venient to be forced to transition to a Np = 5 
scheme because 1) the b-quark clearly plays no 
substantive role in this kinematic range, and 2) 
both the PDFs and ttsifj) will have discontinuities 
at /z = rrif,. 

Clearly it is more reasonable to have the option 
to work consistently in a Np — 4 flavor scheme 
even for fi > mi,. If PDFs were generated such 
that the Np = 4 and Np = 5 schemes co-exist in 
the region /i ^ mi,, then the user could select Np 
by choice. 

The lower portion of Figure [9] illustrates how 
this might be implemented. The PDFs can be 
generated such that the Np scheme is available 
for all /i > mpfp. Thus, for fi — 5 GeV the user 
would have access to schemes with Np = {3,4, 5} 
and can select the scheme by specifying Np in 
addition to {x, Therefore, the user could an- 
alyze their fj, e [2, 5] GeV data set consistently in 
a single Np = 4 scheme, and choose to transition 
to the Np = 5 scheme at a higher n value to be 
specifled by the user. 

The Matching Point: Although the Match- 
ing Point can be set to any ^ value in the region of 
m^p, we shall argue that the choice /zm = mjq^ 
is optimal. 

First, we note that the Matching Point should 
be at or below the Transition Point (/xt > A*m) if 



we desire to avoid downward DGLAP evolution 
(which can be unstable) . Therefore, if we perform 
the matching at the heavy quark mass we have 
the reasonable constraint: > Ma/ = fnNp- 
Second, the matching conditions which deflne 
in terms of /^^ are of the form /^^+^ = 



Aab ^ fb^ with 



A, 



27r 



'Pb- 



In 



Cb~ 



up to 0{al). Here, Pb^a is the DGLAP split- 
ting kernel and Cb^a is a constant]^ The choice 
fiM = misip eliminates the logarithmic terms thus 
simplifying the calculation. 

We also observe that shifting the Matching 
Point from toq to 2toq does not suppress the 
heavy quark PDF as the logarithmic terms com- 
pensate for evolution between mg and 2mQ. 

5. Conclusions: 

The computation of heavy quark production 
has historically been challenging both theoreti- 
cally and experimentally. On the theoretical side, 
the heavy quark introduces an additional mass 
scale which complicates the calculations. On the 
experimental side, the data for heavy quark pro- 
duction has typically differed from the theoretical 
predictions by a signiflcant factor. Recent theo- 
retical developments enable us to incorporate the 
heavy quark mass into the calculation both dy- 
namically and kinematically. These calculations 
have been used to produce matched PDFs in- 
corporating the full mass dependence. Updated 
analyses show improved agreement between data 
and theory for both HERA and Tevatron mea- 
surements. 

Improved experimental precision will demand 
NNLO accuracy from the theoretical calculations, 
and this introduces a number of issues not present 
at the NLO order. There is progress underway on 



"The matching conditions are determined entirely by the 
DGLAP evolution kernels up to a constant term which 
must be computed. At NLO, the constant term is zero 
such that the PDFs are continuous; at NNLO, this term 
is non-zero. 
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both the PDFs and the hard-scattering calcula- 
tions, and this should ensure we are well prepared 
for the upcoming LHC data. 
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